Dynamic scaling in the vicinity of the Luttinger hquid fixed point 
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We calculate the single-particle spectral function A{k,uj) of a one- dimensional Luttinger liquid by 
means of a functional renormalization group (RG) approach. Given an infrared energy cutoff A = 
Aqc"', our approach yields the spectral function in the scaling form, At^{kp -\'P,Lo) = tZiAi{p£^,lut), 
where fc^ is the Fermi momentum, Zi is the wave- function renormalization factor, r = 1/A is the time 
scale and ^ — vp /A is the length scale associated with A. At the Luttinger liquid fixed point {I —f oo) 
our RG result for A{k,u>) exhibits the correct anomalous scaling properties, and for k = zL/cf agrees 
exactly with the well-known bosonization result at weak coupling. Our calculation demonstrates 
that the field rescaling is essential for obtaining the crossover from Fermi liquid behavior to Luttinger 
liquid behavior from a truncation of the hierarchy of exact RG flow equations as the infrared cutoff 
A is reduced. 
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PACS numbers: PACS numbers: 71.10-w, 71.10.Hf 



I. INTRODUCTION 

For many years the renormalization group (RG) has 
been used to study interacting Fermi systems in one spa- 
tial dimension (Id) .El The success of RG methods in Id 
relies on the fact that at low energies, the two-body inter- 
actions between fermions can be parameterized in terms 
of four couphng constants, which are usually called gi 
(backward scattering) , 52 (forward scattering of fermions 
propagating in opposite directions), 173 (Umklapp scatter- 
ing), and 34 (forward scattering of fermions propagating 
in the same direction). Most authors have focused on the 
calculatiaiji.x)f the RG /3-functions, which in the Wilso- 
nian RGclErQ describe the flow of these couplinga^as the 
degrees of freedom are eliminated and rescaled.tl How- 
ever, the RG /?-functions do not completely describe the 
physical behavior of the system. In particular, the RG 
/3-functions do not contain information about the single- 
particle excitations. To investigate these, one should cal- 
culate the momentum- and frequency-dependent single- 
particle spectral function A{k,u!), which at temperature 
T = is related to the imaginary part of the single- 
particle Green function G{k, ui) via 



A{k,uj) = liRG{k,uj + iO) 

TT 



(1.1) 



In the absence of backward and Umklapp scattering, the 
leading asymptotic long-distance and long-time behavior 
of the single-particle Green function G{x, t) in the space- 
time domain can be calculated via bosonization if the 
energy dispersion is linearized arouad. the Fermi points 
(Tomonaga-Luttinger model, TLM)J3'LI Obviously, in or- 
der to obtain the spectral function, one should calculate 
the Fourier-transform G{k,uj) of G{x,t), which in gen- 
eral cannot be done exactly. For the spinless TLM with 
f/2-iiiteractions a mathematically non-rigorous but phys- 
ically reasonable asymptotic analysis yields at tempera- 
ture T = for wave- vectors close to -iikp and for low 



energiesj.""'^""^''^^' 

ylTL(±(fci^ +p),^) 



x\u:-v,pr^+'^'^\^ + VcPy''^ , (1.2) 

where Aq is some ultraviolet cutoff with units of energy 
(for example a band width cutoff). At weak coupling the 
anomalous dimension 77 and the velocity Vc of collective 
charge excitations are to leading order 
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Here vp is the bare Fermi velocity and 
9 = 92/{ttvf) 



(1.3) 
(1.4) 

(1.5) 



is the dimensionless coupling describing forward scatter- 
ing of electrons propagating in different directions. Eq. 
(1.2) satisfies a simple scaling law: for an arbitrary length 
^ we may write 



ATL{±{kF+p),i^) = r[^ 



Ati^{p^,ut) 



(1.6) 



where t = S^/vp is the time scale associated with ^, and 
the length ^0 = vf/^o corresponds to the ultraviolet 
cutoff Aq. The dimensionless scaling function Atl(q, e) 
is 

^tl(<7, e) - |e(e2 - (vqf)\e - vq\-'+^/^\e + vq\^/^ , 

(1.7) 

where q and e are dimensionless variables, and v = Vc/vp 
is the dimensionless velocity renormalization factor. Ob- 
viously, the dynamic exponent (defi ned via r oc ^^) is 
z — 1. The dynamic scaling function (1.7) is scale invari- 
ant, so that for an arbitrary scale factor s 



iTL(sg,se) = s ^+"iTL(g,e) 



(1.8) 
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This scale invariance is a consequence of the fact that the 
TLM represents a critical s*stem, corresponding to the 
Luttinger liquid fixed point Jla Hence, the RG /3- function 
of the TLM vanishes identically. It is generally accepted 
that the TLM describes the generic low-energy and long- 
wavelength properties of one-dimensional Fermi systems 
with dominant forward scattering. Thus, the TLM is an 
effective model which in the regime where backward and 
Umklapp scattering are irrelevant emerges at low energies 
when the high-energy degrees of freedom are integrated 
out in a Wilsonian RG. |— in 

In the theory of dynamic critical phenomenal tI^i it is 
usually assumed that correlation functions close to the 
criti cal point can be written in the dynamic scaling form 
(1.6), with some dynamic scaling function A{p^,ll!t). 
The dynamic scaling function depends on ^ only implic- 
itly via and lot — lo^/vf- For example, at finite tem- 
peratures, where we may identify ^ = vp/T and r w 1 /T, 
the scaling function At-!^{p^,ujt) of the TLM hiia been 
discussed by Orgad, Kivelson, and collaborators^ 

Suppose now that we iterate the RG to some large 
scale factor I — ln(Ao/A), and write the resultin g sp ectral 
function in the scaling form analogous to Eq. (|l.6|). 



Ai{p^,ujt) 



(1.9) 



According to the above argument, we would expect that 
for I ~^ oo the scaling function Ai{q,e) becomes inde- 
pendent of I and approaches Atl(9, e) for large \q\ and 
|e|. In this work we shall show that this expectation is 
not quite correct: even if both |g| and |e| are large com- 
pared with unity, the difference |e| — v\q\ can be small. 
In this regime we find that Aoo{q,£) = limi ^ oo Ai{q, e) 
behaves quite differently from the scaling function of the 
TLM. Physically, this is due to the fact that the spec- 
tral line-shape of a generic one-dimensional Fermi sys- 
tem with dominant forward scattering exhibits in general 
some non-universal features which are determined by ir- 
relevant couplings and thus spoil the universality of the 
spectral function. In fact, even the spectral fipjction of 
the TLM exhibits some non-universal features .llil 

We shall here attempt to calculate the spectral func- 
tion of the f;2-TLM using RG methods. Therefore we 
shall calculate the RG flow of the irreducible two-point 
vertex. Recall that the usual RG /3-function describes 
the flow of the momentum- and frequency-independent 
part of the four-point vertex. Surprisingly, the problem 
of calculating the spectral function of a strongly corre- 
lated fermionic system like the TLM via RG methods has 
not received much attention. Jii fact, with the exception 
of the recent work by Ferrazjlj where the spectral func- 
tion of a special 2d Fermi system has been calculated by 
means of the field theory RG, we are not aware of any 
RG calculation of the spectral line shape of a strongly 
correlated fermionic many-body system. Note that in 
order to obtain a non-trivial k- and w-dependence and 



the anomalous scaling properties given in Eq. (1.6), one 



should retain infinitely many couplings which are irrele- 
vant by naive power counting. In other words, one needs 
to keep track of the RG flow of coupling functions. At 
the first sight, this seems to be a formidable task, which 
is impossible to carry out in practice. Nevertheless, in 
this work we shall show that at weak coupling a simple 
truncation of the hierarchy of functional RG equations 
for the irreducible vertex functions yields an expression 
for the spectral function which has the correct anomalous 
scaling properties and, at least for k = ikp, agrees with 
the exact result known from bosonization. 



II. EXACT FLOW EQUATIONS 

Originally, exact RG fiow equations have been devel- 
oped in field theory and statistical phys ics , t o ^tudy . sys- 
tems in the vicinity of a critical point .llIHl30'E2lEil Re- 
cently, several authors have started to apply thesa^ieth^ 
ods to interacting Fermi systems at finite densities .E3lH3'0 
In the normal state the existence of a Fermi surface adds 
some new complexity to the problem: because the single 
particle Green function exhibits singularities on the en- 
tire Fermi surface, the Fermi surface plays the role of a 
continuum (in c? > 1) of RG fixed points. However, this 
fixed point manifold is not known a priori, and should be 
calculated self-consistently within the RG. In Ref. ^ we 
have shown how the true Fermi surface of the interacting 
system can be calculated from the requirement that the 
RG approaches a fixed point. Fortunately, in Id we know 
a priori that the Fermi surface (which consists of two. 
points ikp) is not renormalized by the interactions ,Ej 
so that the above mentioned problem of a self-consistent 
determination of the Fermi surface does not arise. 

Even if the shape of the Fermi surface is known, the 
implementation of the exact RG for Fermi systems is 
non-trivial and some open problems remain. In par- 
ticular, in Ref. we have pointed out that the RG 
approach proposed in Refs. p2| , p3|p^ includes only the 
mode-elimination step, but completely omits the usual 
rescaling. We have further speculated that this incom- 
plete RG transformation might be responsible for the 
runaway flow to strong coupling which is typically en- 
countered in these works. The purpose of this article is 
to further substantiate this claim by applying the func- 
tional RG method to a non-trivial yet exactly solvable 
many-body problem where the results can be critically 
compared with known results obtained by other meth- 
ods. 



A. Definition of the model and notations 

Starting point of our calculation are the exact fiow 
equations derived in Ref. ^ which give the RG fiow of 
the one-particle irreducible vertices when degrees of free- 
dom are integrated out and a subsequent rescaling step 
is applied. For convenience we shall now briefiy present 
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these equations for the special case of Id, using a shghtly 
different notation than in Ref. We shall work at zero 
temperature, but expect that our results also describe 
temperatures T > as long as we consider energy scales 
large compared with T. In Id, the Fermi surface con- 
sists of two points, akp, where a — ±1 labels the right 
and left Fermi point. An arbitrary wave- vector k can be 
written as 



k = a{kF + p) 



(2.1) 



where a = -1-1 for fc > and a = — 1 for fc < 0. Note that 
the deviation p = a{k — akp) from ifci? is always mea- 
sured locally outwards, which is differeijJ_^aiii the con- 
vention usually adopted in bosonization.ErQu'EJ For small 
p we may expand the non-interacting energy dispersion 
Cfe around the Fermi points. 



P 

<^a{kF+p) = '^'kp + VfP +7, ^ 



(2.2) 



Here m is the bare mas s of the fermions. As discussed in 
Ref. I25I, the expansion (2^) is around the true Fermi mo- 



mentum of the interacting many-body system. Of course, 
in \d the interacting and non-interacting kp are identi- 
cal if we compare systems with the same density. Note 
that by time reversal invariance ej, — e_/j, so that vp and 
771 are independent of a. For \p\ <C kp it is reasonable 
to linearize the ene rgy dispersion and neglect the term 
of order p^ in Eq. (2.2). Introducing an infrared cutoff 
A and an ultraviolet cutoff Aq with units of energy, the 
non-interacting Matsubara Green function is 

r>o ( (u ^ \ ■ \ Q(Ao > vp\p\ > A) 
GA,Ao(a(^F +p),«t^) = : , 2.3) 

' " tUJ — Vpp 

where ico is a fermionic Matsubara frequency, and 



Q(X2 > X > Xl) 



Later we shall also write 



8(x2 — x) — Q{xi — x) 

1 if X2 > X > Xl 

else 



(2.4) 



&{X2 > Xl) = Q{X2 -xi)^{l'^ ""l^l^ ^1 . (2.5) 



Taking the limit Aq — > cxo in Eq. (2.3) amounts to ex- 
tending the linear energy dispersion on both branches 
from —00 to -|-oo. The unphysical states with energies 
far below the Fermi energy introduced in this way are 
occupied according to the Pauli principle. It is generally 
accepted that this filled Dirac sea is dynamically irrele- 
vant and does not modify the low-energy physics. How- 
ever, the precise way in which the cutoff is removed can 
affect the nvKiierical value of the various Luttinger liquid 
parametersO At this point we shall keep Aq finite and 
assume that the bare electron-electron interaction of the 
model is characterized by a momentum- and frequency- 
independent totally antisymmetric irreducible four-point 
vertex of the form 



^^Koi^'l^ ^2, Ki) = ^aia^;a2ai50 , 



(2.6) 



where K = (A:, iuj) = (a{kF + p)-, ^^), and 



r) 771 



(2.7) 



is antisymmetric with respect to the exchange ai <-> q;2 
or a'^ <-> a'2- For later convenience we have introduced 
the notations 



^a'a';a2ai = ^a' a^Sa' a^, (direct term) (2.. 



(exchange term) . (2.9) 



In the limits Aq — + oo and A the Green function of 



the model defined by Eqs. (2.3) and (2.6) can be calcu- 
lated exactly in the space-time domain via bosonization 
(TLM with g2-interactions). The generally accepted re- 
sult for the spectral function is given in Eq. (1.2)Bj 



B. Scaling functions and flow equations 

Suppose now that we reduce the infrared cutoff by set- 
ting 



A = Age" 



(2.10) 



and follow the evolution of the correlation functions as 
the logarithmic flow parameter I increases. Noting that 
the infrared cutoff A has units of energy, we may define 
a length scale ^ and an associated time scale r. 



i/A = e/«F 



(2.11) 



From these quantities, we may construct dimensionless 
scaling variables, 

q = p£, = (ak-kp)^ , e = cjt = cj/A , (2.12) 

and write the exact single-particle Matsubara Green 
function of the theory with infrared cutoff -A-ip-Rd ultra- 
violet cutoff Aq in the dynamic scaling forinlljo 



GA,Ao{k,iuj) = tZiGi {p^,iu!T) 



(2.13) 



where Gi{q,ie) is a dimensionless dynamic scaling func- 
tion. Here the wave-function renormalization factor Z; is 
related to the irreducible self-energy Y,\{k,iijj) as usual. 



Zi = 



1 



1 



d{iuj) 



(2.14) 



j=0 



Note that by time-reversal symmetry Zi and Sa are in- 
dependent of a. It is also convenient to set 



Gi{q,ie) 



e(e' > |g| > 1) 



where 



ri{q,ie) = Zi[ie - q] {q,ie) 



(2.15) 



(2.16) 
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and 



We have introduced the notations 



(2.17) 

Here iw) = HniA— >o 5]A(fc, iu) is the exact self-energy 
of the model without infrared cutoff. Following Ref. |2^, 
we also define the scaling functions for the higher order 
irreducible vertices, 



f I \Qi) ■ ■ ■ 1 Qn'j Qn 



i) 



v^^-'K^^-^ZrTf'\K[, ... ,K'„;K„,... ,K^), (2.18) 



where T 



A 



K' ■ K 



, Ki) are the usual one- 



particle irreducible 2n-point vertices, vt^ = {■KVp)~^ is 
the density of states of non-interacting spinless electrons 
in Id, and K = [k^ioj) and Q = {a,q,ie) are composite 
labels. 

T he re scaled irreducible two-point vertex defined-Jn 
Eq. (2.17) satisfies the following exact flow equation,Ea 



dir[^'{Q) = {i-m-Q-dQ)r\^>{Q) + r\''{Q) , (2.19) 



(2), 



where 



fp)(g) 



de' 



2 



X r!"(<3,0';<3'.<3), 

and r]i is the flowing anomalous dimension, 

d,Zi 



-dilnZi 



'A 



(2.20) 



(2.21) 



Q ■ dg ^ qdg + ede , 



Gi{Q) = 



n{Q) 



(2.22) 



(2.23) 



Eq. ( ^.19 ) is equivalent with the following integral equa- 
tion. 



I 



+ / die*-'''(*)(l-r;i_Of[iUe"*Q) 



+ / die*-'''(*)r|!.^i(e~*Q) 



(2.24) 



where e *Q = (a, e *q, e *ie), and we have defined 



ili{t) = / dt'rii^t+t' 
Jo 



The flow of the in homogeneity (Q) on the right-hand 
side of Eq. (2.19) is controlled by the scaling function of 

the four-point vertex f[^\ which satisfies the exact flow 
equation 



(2.25) 



dir\^\Q[,Q'2;Q2,Qi) 



r\^\Q[,Q2;Q2,Qi) + t[^\Q'i,Q2;Q2,Qi) , (2.26) 



where 



t['\Q'i,Q2-:Q2,Qi) = -IY. J'^iJ ^GiiQ)r['\Q[,Q'2,Q;Q,Q2,Qi) 



1 



de 



'^^y ^ [Gi{Q)Gi{Q') + Gi{Q)Gi{Q') 

[f (Q'l , Q2;Q', QWl'^ {Q,Q';Q2,Qi) 



1 
2 

f (Q'l, Q'; Q, Qi)f f ^ (^2, Q; Q', Q2) 



r[^\Q2,Q':Q,Qi)^[^\Q'i,Q;Q',Q2) 



K'=Kt+K2-K 
K'=K+Ki-K[ 
K'=K+Ki-K' 



(2.27) 



Here vf^ is the dimensionless irreducible six-point vertex. As Eq. ( p.2C ) can be converted into Eq. (2.24), Eq. (2.2i) 
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can be transformed into the integral equation 



f|'^(Q'i,Q2;Q2,Qi) = f|foWi'Q2;Q2,Qi)- f dte-'^^^*hr,i_tr\%{e~'Q[,e''Q'^;e-'Q2,e''Qi) 

Jo 



+ / dte 
Jo 



(2.28) 



C. Classification of couplings in the hydrodynamic 
regime 

The above functional RG equations determine the RG 
flow of the momentum- and frequency-dependent irre- 



h). By ex- 



ducible vertices r\'^"\Q[, ... ,Q'n]Q 
panding these vertices in powers of the dimensionless 
scaling variables qi = piS^ and = w^r, we obtain the 
RG flow of the coupling constants of the model, which 
are the coefficients in this multi-dimensional Taylor ex- 
pansion. Following the usual terminology, couplings with 
positive scaling dimension are called relevant, couplings 
with vanishing scaling dimension are called marginal, and 
couplings with negative scaling dimension are called irrel- 
evant. The irrelevant couplings grow at short distances 
(i.e. in the ultraviolet) and spoil the renormalizabiliy of 
the theory. In this case the usual field theory RG is not 
applicable. On the other hand, in a Wilsonian RG the 
irrelevant couplings can be treated on the same footing 
as the relevant ones. For this reason the Wilsonian RG 
is very popular in condensed matter physics, where one 
usually has a physical ultraviolet cutoff and renormaliz- 
ability is not necessary. 

What is the regime of validity of the above expansion? 
Obviously, the expansion of the rescaled vertex functions 
to some finite order in powers of the dimensionless scaling 
variables qi = pi^ and = uJiT can be justified if these 
variables are small compared with unity. As long as the 
infrared cutoff A = vp/^ — 1/t is finite, this condition 
can be satisfied for sufficiently small wave-vectors and 
frequencies. If we identify ^ with the correlation length 
and T with the corresponding characteristic time scale, 
then the expansion of the f ^"^(Q^, ... ,Q'„;Qn, . . . ,Qi) 
in powers of the scaling variables qi = and = w^r 
can be viewed as an expansion of the dynamic scaling 
functions in the hydrodynamic regime^ where one is inter- 
ested in length scales larger than ^ and time scales longer 
than T. Of course, a normal metal at zero temperature 
is a critical system, where ^ and t diverge. Close to this 
critical state (where ^ and r are large but finite), we ex- 
pect that the correlation functions assume some scaling 
form. This corresponds to the critical regime (also called 
scaling regime). For example, we may describe a normal 
metal at finite temperature in terms of flow equations 
derived for T = if we stop scaling the flow equations at 
a finite length scale of the order of ^ « vp/T. 

For certain systems it may happen that the couplings 
defined in the hydrodynamic regime smoothly evolve into 



analogous couplings in the scaling regime, where both 
qi and ct are largeu-This is the essence of the dynamic 
scaling hypothesis .t3 For example, for a Fermi liquid we 
expect that such a smooth crossover between the hydro- 
dynamic and the scaling regime indeed exists, because 
by definition the self-energy is analytic and hence can be 
expanded in powers of momenta and frequencies. On the 
other hand, the two-point function of a Luttinger liquid 
is known to exhibit algebraic singularities, so that there 
cannot be a smooth connection between the hydrody- 
namic and the critical regime. In this case one expects 
the analytic properties of the dynamic scaling function 
F; {q, ie) to change as we move from the hydrodynamic 
into the scaling regime. In this work we shall show that 
this is indeed the case and give an approximate expres- 
sion for the scaling function. 

Let us now classify the couplings according to their 
relevance in the hydrodynamic regime. First of all, the 
momentum- and frequency-independent part of the two- 
point vertex. 



(2.29) 



is relevant. Her e Fp-* (0) stands for r\^' (a, q — 0,ie = iO) 



(2), 



From Eq. ( p.l9D it is easy to see that /i/ satisfies the exact 
flow equation 



difii = (i-?//)A/ + rp^(o) . 



(2.30) 



Note that in general we have to fine-tune the other cou- 
plings such that the relevant coupling jli approaches a 
fixed point. The fixed point value of fli for Z — > cxd is then 
determined by 



(i-??oo)/ioo = -r£'(o) 



(2.31) 



Because Fjxj (0) implicitly depends on jlac, Eq. (2.31) is 
really a self-consistency equation for the fixed point value 

of ill . 

There are two marginal couplings associated with the 
irreducible two-point vertex F, (Q): the wave-function 
renormalization factor Zi , and the Fermi velocity renor- 
malization factor w/. Using the definition ( 2.14 ), it is easy 
to show that Zi can be expressed in terms of our rescaled 
two-point vertex as follows. 



Z, 



dt\^\a,0,ie) 



d{ie) 



(2.32) 



e=0 
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The dimensionless Fermi velocity renormalization factor 
can be written as 



Z, 



(2.33) 



q=0 



Note that for a Fermi hquid vi can be identified with the 
dimensionless inverse effective mass renormalization, 



m 
ml 



(2.34) 



where m is the bare mass and mj" is the effective mass of 
the model with infrared cutoff Aoe~'. The expansion of 
the scaling function associated with the two-point vertex 
for small q and e then reads 

ff\Q) = fii + (l - Zi)ie + [Zi - vi)q 

+ 0{q\e\qe), (2.35) 

so that for small q an d e th e dimensionless inverse prop- 
agator defined in Eq. ( ^.16 ) is given by 



n{Q) = Zi[ie-q]+Tf\Q) 

= ze - viq + fH + 0(<Z^ e2, qe) . (2.36) 

The exact flow equations for the mar ginal co uplin gs Zi 
and Vl are easily obtained from Eqs. ( ^.19 ), (2.21) and 

(Pi, 



diZi 



-mZi , 



dt[^\a,q,tO) 



dm 



dq 



(2.37) 



(2.38) 



9=0 



Differentiating Eq. (2.19| ) with respect to ie, setting e = 0, 
and using Eq. (2.3S), we see that the flowing anomalous 

* (2) 

dimension rji is related to the function Fj (Q) via 



(2.39) 



€=0 



* (2) 

Since Fj {a,0,ie) is defined in terms of the four-point 



vertex f[*^ via Eq. (|^), we can obtain the following 
explicit expression for the flowing anomalous dimension, 



m = 



2 ^ 



dq' 



2^ 



GiiQ') 



dr['^\a,0,ie,Q';Q',a,0,ie) 



d{ie) 



e=0 



This expression is also valid in dimensions d > 1 pro- 
vided the discrete index a is replaced by a d-dimensional 
unit vector which labels the points on the Fermi surface. 



Thus, for a given irreducible four-point vertex, Eq. (2.4C) 



can be used to directly calculate the anomalous dimen- 
sion of any normal Fermi system. 

There is one more marginal coupling, namely the irre- 
ducible four-point vertex at vanishing momenta and fre- 
quencies. Because by construction f ["'"'(Qi, Qa, Qi) 
is antisymmetric with respect to the permutation of the 
incoming or the outgoing ferniions, for vanishing external 
momenta and frequencies the four-point vertex must be 
of the form 

f|''^(ai,0,a2,0;a2,0,ai,0) = Aa'^a'^-a^aiQi , (2-41) 



where ^05^0^:0201 is defined in Eq. (2.7). Hence, the 
marginal part of the four-point vertex can be expressed 
in terms of a single marginal coupling 5;, which satisfies 
the exact flow equation 



digi = -2r]igi + Bi . 
Here Bi is defined by 

rr^(a'i,0, Q;2,0;a2,0,ai,0) = A^'^a'^-o 



,Bi 



(2.42) 



(2.43) 



with f f ^ given in Eq. ( ^^ ) 



It is convenient to subtract the relevant and marginal 
couplings from the exact vertices, thus defining coupling 
functions containing only irrelevant couplings. The irrel- 
evant part of the two-point vertex is 



^{2-^l,z,v) 



(Q) = ff\Q)~TY>(a,0,iO) 



I 



die 



dff\a,q,iQ) 



9=0 



fii - ie{l - Zi) - q{Zi - Vl) , (2.44) 



where the superscripts indicate the marginal and relevant 
couplings to be subtracted. T he re scaled inverse propa- 
gator ri[q,ie) defined in Eq. (2.16) can then be written 
as 



ri{q, ie) = ie - viq + /i( + f 



{2-fj.,Z,v) 
I 



(Q) . (2.45) 



Finally, the irrelevant part of the four-point vertex is 



=^(4-3) 
I 



(Q'l, Q2; Q2, Qi) 



f J'^(0'i, Q2; Q2, Qi) - f r(Q'i, Q^; Q2, Qi) 

f| ^(Qi, Q2; Q2, Ql) ~ ^Q'jQ^;a2Ql5i • 



(4)/ 



9i=e;=0 

(2.46) 



The exact flow equations for the subtracted vertices are 
formally identical with the flow equations fo r th e un- 
subtracted vertices given in Eqs. ( I2I9I ) and ( |2^ ), ex- 



fay 

cept for the fact that the inhomogeneities Fj (Q) and 
(2.40) ^i^\QiiQ'2'^Q2,Qi) on the right-hand sides should be 



replaced by their subtracted versions. 



(Q) = fp)(Q)-Fj^^(a,0,zO) 
dr[^\a,q,iO) 



f(2). 



I 



afp^(a,0,ie) 



die 



c=0 



dq 



(2.47) 



9=0 



t'i^~'\Q'i,Q2;Q2,Qi) 

f(4) 



t[^\Q[,Q2;Q2,Qi) 



-T\''>iQ[,Q'^;Q2,Qi) 



(2.48) 



III. APPROXIMATE SOLUTION OF THE 
FUNCTIONAL FLOW EQUATIONS FOR THE 
TWO-POINT VERTEX 

So far no approximation has been made, except for the 
linearization of the energy dispersion. In order to make 
progress, we need to truncate the hierarchy of flow equa- 
tions. We now give a truncation scheme which at weak 
couphng reproduces the known scahng properties of the 
spectral function A{zt{kp +p),uj) of the TLM. Moreover, 
at least for p = 0, we recover the exact weak coupling re- 
sult for the spectral function known from bosonization. 



A. One-loop flow 

To begin with, let us briefly discuss the flow of the 
relevant and the marginal couplings within the one-loop 
approximation. We first consider the t wo- po int vertex. 
The exact flow equation is given in Eq. ( 2.19 ). The flow 
is coupled to the irreducible four-point vertex via the 
inhomogeneity (Q) on the right-hand side, as given 
in Eq. ( ^.20 ). Within the one- loop approximation, it is 
sufficie nt to approximate the four-point vertex appearing 
in Eq. (E2O) by 



t'i^\Q,Q';Q',Q) = 6^^^^,gi + 0{~gf) 

where we have used the fact that Aaa'-a'a = 
Sa,^a'- In this approximation 



(3.1) 



-| [Qifii + ii) + e(A/ - ii)] + 0{gf ) . (3.2) 



For 1/2; I < I'D; I this simplifies to 



fp)(Q) 



m 
2 



Oi~9f 



(3.3) 



As shown below, in the physically relevant case ui = 
Ojgi ), wh il e vi is of the order of unity. From Eqs. ( 2.30| ) 
and (2.37)-(2.39) wc conclude that for small gi 



dm 
diZi 

dm 



M; - 

0, 
. 



m 
2 



(3.4) 

(3.5) 
(3.6) 



Finally, let us consider the momentum- and frequency- 
independent part of the four-point vertex. Using the fact 



that r[^^ is of order gf an d hence can be neglected 
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we obtain from Eq. (2.27) to second order in gi 

a 

^'12" ^"^^'^ -.a.' a-^aa' ;a2Q:i 

X {Gi{Q)GiiQ') + Gi{Q)Gi{Q')j,j,,^K^^j,^_j, 

X {GiiQ)GiiQ') + GiiQ)Gi{Q'))]^,^^^j^^_j^,} . 

(3.7) 

Because on the critical manifold fli = 0(gi). we may 
approximate on the right-hand side of Eq. (3.7) 



le - viq 



le - viq 



(3.8) 



The integrations in Eq. (3.7) can then be performed. We 
obtain for the final result 



-~9f 



Xi{qi - q2,iei + K2) 
+ -a2X/(9i +9i,*ei -ie'i)] 



Here 



X;(g,*e) = e(e'-1> |<7|) 



ri,(g,ie) = e(e' + l>|g|>2) 



^K2 + |g|) 

£2+^2(2+1^1)2 
So 



le + viq 

where we have introduced the notation 
Sg sign(g) . 



(3.9) 

(3.10) 
(3.11) 

(3.12) 



Setting external mo men ta and frequencies equal to zero 
we obtain from Eq. (13. 



f|'^^(ai,0,Q;2,0;Q;2,0, ai,0) = -gfAa'^a'^.a2ai 

x{(5a,,„,X!(0,jO)-5a,,„,ri,(0,iO)} = 0. (3.13) 

Hence, the coefficient of order gf in the wea k cou pling 
expansion of the function Bi defined in Eq. ( 2.43 ) van- 
ishes. Anticipating that rji — 0{gf), we obtain from Eq. 
( ^.42 ) within the one-loop approximation 

di~gi = . (3.14) 
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In this case jli = ^ for all I, so that we obtain a RG 
fixed point. The one- loop flow equations given above 
have also been discussed by Shankar,L2l who derived these 
equations within the framework of the conventional mo- 
mentum shell technique. Because within the one-loop ap- 
proximation Z; = 1, this approximation is not sufficient 
to detect the non-Fermi liquid behavior of our model. 



FIG. 1: O ne-loop RG flow in the g — /i-plane, see Eqs. (3.4) 
and (3.14). The thick solid line is a line of fixed points. 



Two-loop flow of the two-point vertex 



The vanishing of the /3-function of the marginal coupling 
gi of the TLM is well-known. However, the one-loop RG 
fiow of the relevant coupling //; is non-trivial, because 
the initial value of fli has to be fine-tuned such that fli 
does not exhibit a runaway fiow. As emphasized in Ref. 

the requirement that fli fiows into a fixed point is 
equivalent with the statement that the initial kp is the 
correct Fermi mcirnentum, in agreement with the Lut- 
tinger th eore m.P^P ^ The flow in the g — /i-plane implied 
by Eqs. (3.4) and (3.14) is shown in Fig. ^ Obviously, 
at the one-loop level the marginal couplings gi, Zi and 
vi do not flow, while the relevant coupling fli exhibits a 
runaway flow unless the initial value is chosen such that 



90 
2 



(3.15) 



It is now straightforward to calculate the full two-point 
vertex within the two-loop approximation and to study 
the emergence of Luttinger liquid behavior when the in- 
frared cutoff is reduced. The crucial observation is that 
for large values of the flow parameter I the momentum- 
and frequency-dependent part of the four-point vertex 

f P ^\Qi, Q'2:Q2,Qi) can be expanded in powers of the 
marginal coupling gi. This follows from the fact that 
the RG trajectory approaches the manifold defined by 
the relevant and marginal couplings, so that the irrele- 
vant couplings, bec ome local functions of the relevant and 
marginal onesJl3c3 

To calculate the two-point vertex at the two-loop or- 
der, we need to know the momentum- and frequency- 
dependent part of the four-point vertex within the one- 
loop approximation. By definition, 



r[''^(Qi, Q'2', Q2, Ql) — Aa'^a'2;a2ai9l + f J" \Q'i, Q'2'1 Q2, Ql) 



(3.16) 



where to order gf we obtain from Eqs. (2.28) and (3.9), 



Here 



fp"^^(Ql,Q2;Q2,Ql) ~ f|io''^(Qi,Q2;Q2,Ql) - ~gi\Aa\a'^,a.2a^xM\ - 92,»ei +»e2) 
+ Ea'^a'^-a-ioc^ ^a^^Ml - q'i,i(^l ^ ie'i) + 5at.-a2Xl{qi + q'l,i<^l - ^e'l)] 



(3.17) 



H,(g,ie) = n,,o(g,ie) , «e) = X;,o('?, ie) , 

where for arbitrary initial value ^ 0; 

X/Jo(9:«e) = / [x;-t(e^*g,e~*ie) - x/_((0,zO)] 



-^|(?-Me(|g|>e'-e'°)-f iln 



vf{2 + e-~''°\q\Y + 
«2(2-e-'|9|)2 + e-2'e2 



^2 I g^2io ^2 



e(e' - e'" > M) j> , (3.19) 



(3.18) 
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and 



n;,Zo((?,ie) 



dt 



1 Sq 

2 ie + viq 



q\ ~ 2e^°)e{e^ + e'° > \q\ > 2e^°) + (2e' - k|)e(2e' > \q\ > e' + e'") >> . (3.20) 



Note that Eq. (3.17) is manifestly antisymmetric with 
respect to the exchange of the incoming or the outgoing 
labels. For the calculation of the two-point vertex we 
only need the irreducible four-point vertex at Qi = Qi 
and Qi = Q'^, i.e. 



(PB) channels, we have 

(3.22) 



f[^)(Q,Q';Q',Q) = 5,,_„,5j + 5,,„,5fnj(g-g',ze-ze') 



Xi{q + q',ie - ie') + xM - q',ie + ie') 



(3.21) 



A graphical representation of this equation is shown in 
Fig.|. We recognize the usual perturbative contributions 
to the effective interaction: the first term of order on 
the right-hand side of Eq. ( ^.21 ) corresponds to the con- 
tribution from the zero-sound (ZS) channel, the second 
term corresponds to the Peierls channel (sometimes also 
denoted by ZS' channel), and the last term corresponds 
to the BCS channel. Note that the zero-sound channel 
gives rise to an effective retarded forward scattering inter- 
action of the (74-type (using the g-ology language) , which 
is generated by integrating out degrees of freedom. 

Given an approximate expression for the four-point 
vertex, we may calculate the inhomogeneity (Q) on 
the right-hand si de of the fl ow eq uation for the two-point 
vertex, see Eqs. (2.19) and ( ^.2C| ). Separating the contri- 
butions due to the zero sound (ZS) and the Peierls-BCS 



where 

p(2,ZS) 



(<?, «e) 

^Q{e^>\q\>l). ''^'"^ 



vi{2 + \q\) + iesq 

e(2e'-l>M>e')-|^^ld^ 
vi(2 + \q\) + leSq 



(3.23) 



p(2.PB 



9f 



I (g,ie) = -y / dq 



, f de'5{\q'\-l) 



~2 

9i 



27r ie' — viq' 
X [xi{q + q',ie-ie') - ^((q - q', + ie')] 

vi{A - \q\) + ieSq 



^L.Sq\Qil> |g| >2-e')ln 



w/(2e' + \q\)+iesq 



e(e' > \q\ > l)ln 
e(e' - 2 > |g|)ln 



vi{2+\q\)+i€ 



viC^e' + 2 - \q\) + ie 
Vi{4:+ \q\) - ieSq 
w;(2e' - \q\) - ieSq 



(3.24) 



Expanding to first order in q and e, we obtain 



+ 




+ 



+ 



-< 

ZS 

< I < 1 < 



BCS 



Ff'™)(g,ze) = g2 0(e'-2) (1 - 2e-')*e 



5? 



-{l + 2e-')viq + 0{e\q',eq) .(3.25) 



Because vf'''^^\q,ie) vanishes for \q\ < 1, it does not 
contribute to th e flow of the marginal couplings Zi and vi . 
From Eq. ( 2.39 ) we obtain for the anomalous dimension 
at scale I, 



m = ^eie'-2){l-2e-'), 



(3.26) 



FIG. 2: Feynman diagrams contributing to the efFecti ye in - 
teraction within the one-loop approximation, see Eq. ( 3.21 ). 
The solid (dashed) arrows represent the fermion propagators 
with momenta close to the left (right) Fermi point, and the 
wavy lines represent the flowing coupling constant gi . 



and from Eq. ( ^!38| ) we find for the flow of the Fermi 
velocity renormalization factor 



dm = |-e(e' - 2)e-' 
2vi 



(3.27) 
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; - 



0.5 











= 0.1 


/ 9o 


= 0.5 



FIG. 3: Anomalous dimension 77; as a functi on of the loga- 
rithmic flow parameter see Eqs. (|3.26|) and (3.27). 



Using the fact th at qi = g \s independent of the flow 
parameter Z, Eq. ( 3.27 ) can be easily integrated, 



.9' 



l + e(e'-2)i^(l-2e-') 



1/2 



Note that vi rapidly approaches a constant Vo 
I, which is to leading order in g given by 



va + -^ + 0{g 



(3.28) 



for large 



(3.29) 



In Fig. y we show the flow of the anomalous dimension 



rji given in Eq. (3.26) as a function of the logarithmic 
flow parameter I. Obviously, rji vanishes for I < In 2, and 
approaches a constant rj for large I, which is given by 



V - 



8vo 



(3.30) 



The above fixed point value of i] agrees with the weak 
coupling expansion of the bosonization result for the 
anomalous dimension of the spinless 52-TLM. Note, how- 
ever, that the exact solubility of the TLM relies on the 
filled Dirac sea associated with the linearized-finergy dis- 
persion. As discussed by Schulz and ShastryjHZl the intro- 
duction of the Dirac sea leads to finite renormalizations 
of the fixed point values of the Luttinger liquid parame- 
ters (such as the anomalous dimension), so that the value 
of 77 beyond the leading order g is modified. This can be 
explicitly verified within our RG approach, where the in- 
troduction of the Dirac sea corresponds to the removal of 
the ultraviolet c utoff Ap. T hen w e should take the limit 
" ^ 00 m Eqs. ( |3.23| ) and ( [3. 24] ), so that the exphcit l- 



dependence on the right-hand sides of these expressions 
disappears. In this case the marginal coupling vi is not 
renormalized at all (divi = 0), so that Voo — vo- Fur- 
thermore, the initial flow of the anomalous dimension iji 
discussed above is also removed by the Dirac sea, so that 
r]i is eff ective ly replaced by its asymptotic limit 77 given 
in Eq. (3.30). To simplify the following calculation of 
the spectral function, we shall from now on work with a 



Dirac sea, which amounts to t aking the limits Aq ^ 00 
and Z — > 00 in Eqs. (|^ and (U^. Then we obtain 



f(^-^s)(g,ze) = -e(M > 1)2,7.,— -^L^— 

uo(2 -I- \q\) + lesq 



Q{\q\ > 1)775, 



3vo\q\ + ieSq 



vo{2 + l^l) + iesq 



, (3.31) 



r(^'P«)((7, te) = 2775, |e(|g| > 1) ln[^o(2 + \q\) + teSq] 
+ 6(1 > |(7|)ln[^o(4-|<7|)+ies,] 
-ln[So(4+k|)-ies,]| . (3.32) 

The correspo nding subtracted function F^, ^'^'^\q) de- 
fined in Eq. (2.47) is in this approximation given by 



'^^^'^''^H'?,*^) = 77s,|2e(|<7| > l)ln[«o(2+|g|) + zes,] 



29(1 > |g|)ln[^o(4-|g| 
2 In [So(4+ \q\)~i€Sq] 
ivn\q\ 



+ e(|g| > 1) 



1 - - 



^o(2+k|) 



- ieSq 
(3.33) 



So far we have calculated the function f ^ ^'^'"''(Q) de- 



fined in Eq. (2.47) perturbatively to second order in g. 
Let us for the moment proceed further within pertur- 
bation theory. In this case we may set 77; = on the 
right-hand side of the exact flow equation ( ^.1S| ) for the 
two-point vertex, which amounts to setting fji{t) — in 



Eq. ( |2.24| ). Assuming that initially r|^o'''^'"^ (Q) = 0, we 



obtain in this approximation 



:,(2-Ai,Z,u) 



- I 



{q,ie) « dte*r 



(e *g, e *7e) 



(Ag, Me) 



(3.34) 



Subst ituting the perturbative result ( |3.33 ) into Eq. 
( |3.34 ), going back to unrescaled variables p — Aq/vp, 



u = Ae, Fp^(Q) = -f- [Ea(/c, itu) - j:{akF,iO)] (see Eq. 
( ^.17| )), and finally taking the scaling limit I — > 00 (i.e. 
A — > 0), we recover the well-known perturbative self- 
energy of the spinless g2-TLM, 



I](a/cF + cep, ioj) — Ti(akF, iw) 



VFP+ 777 - vfp) In 
Id 



It is also easy to check that in the regime 
K |g| < e' , 1 <C |e| < e' , 



we may approximate 



-p,(2-^,Z,t)) 



(g, ie) « 



(2-/^,Z,t)) 
/ — ^00 



(9, «e) 



,(3.35) 



(3.36) 



(3.37) 
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so that the only /-dependence of the self-energy enters via 
the scaling variables q — p£, = P^oc' and e = wr = oj/A — 
we'/Ao- Thus, in the scaling limit the spectral function 
can be written in terms of a scaling function which does 
not explicitly depend on the logarithmic scale factor I. 

We now propose a simple procedure to go beyond per- 
turbation theory: from Eq. (2.24) we know that the exact 
two-point vertex for I 3> 1 satisfies 



t[^-^-'^-''\q,ie) = / dte(i-'')*r(^~'^-^^^)(e~*g,e-**e) 
= [ (iAA-2+')fg-^'^'")(Ag,Aie) , (3.38) 



which differs from the perturbative expression ( [3.34 ) be- 
cause the fixed-point value of the anomalous dimension 
appears on the right-hand side. We now approximate 
the flow function t'^~'^''^'^\q,ie) on the right-hand side 



of Eq. (3^) by 



rg-^-^^^'\q,ie) « 7^s,}^2ei\q\>l)\n[v{2+\q\) + ^es,] 

+ 2e(l> |g|)ln[^(4-|g|)+zes,] 
-2 In [viA+\q\)-ieSq] 

-9(M>i) -,fi^i';!^.- 

v{2 + \q\) + lesq 

- {iesg - v\q\)Y (3.39) 

This expression is alm ost identical with the perturbative 
two- loop result ( 3.33 ), except that the term 775^8 (|g| > 1) 
has been omitted, and the factor vq has been replaced by 
the perturbative Fermi velocity renormalization factor 



~ 1 9 1 

V = I = 1 — 77 



(3.40) 



see Eq. ( |l.4| ). Note that with in p erturbation theory 
the term r]Sqld{\q\ > 1) in Eq. (3.33) is responsible for 
the finite Fermi velocity renormalization given by the 
first term —^vpp on the right-h and side of Eq. (3.35). 
By substituting vq ^ v in Eq. ( [3.33| ), we have implic- 
itly replaced the bare propagator in the scaling limit by 
[ie—vq]~^, thus taking the Fermi velocity renormalization 
in the scaling limit (where the dimensionless variables p 
and e are large compared with unity) self-consistently 
into account. Note that according to Eq. (2.34) the di- 
mensionless quantity v can be identified with the renor- 
malization of the inverse effective mass, which remains 
finite in Id. 



Given Eqs. ( 3.38 ) and ( 3.39 ), we may calculate the 
complete dynamic scaling function for the irreducible 
self-energy, describing the change in scaling behavior as 
we move from the hydrodynamic into the scaling regime. 
In the hydrodynamic regime the scaling function is ana- 
lytic in q and e. On the other hand, for large \q\ and |e| 
the scaling function exhibits algebraic singularities cor- 
responding to Luttinger liquid behavior. In the scaling 



regime (3.36) we i nay use the app roximation (3.37) and 
obtain from Eqs. ( |3.38| ) and ( 0^ ) 



X / 



leSq - v\q\ 



f - 



iesq + v\q\ 



2v 



J I _ij£q_^ v\q\ 



Ai\q\ 



-f[- 



ieSq + v\q\ 



2i\q\ 



(3.41) 



where the complex function f(z) is defined by 
A" 



Z2i^i(l,l + ?7,2-H?7;z) . 

(3.42) 



Here, 2Fi (a, h , c; z) is a hypergeometric function. Note 
that Eq. ( 3.41 ) is restricted to the scaling regime defined 
by ( 3.36| ). To obtain the spectral function, we analyti- 
cally continue f ^ ^'^'"-^ (g, ie) according to ie ^ e + iO. 
Using 



Im/(a; ± iO) = ±7re(x > l)a;-'' , 

and assuming that — u^q^j ;2> 1, we obtain 

Inif(^-^'^'"'(g,e + iO) = 27r77 

€ — vq 



(3.43) 



x<^ e(esq > v\q\) 



4w 



-I- e(-u|g| > esg > -~'iv\q\) 



Q{-?.v\q\ > eSq 



e — vq 



Av 



e + vq 



2w 







e — vq 


Av 





(3.44) 



Mathematically the imaginary part of f m"^''^'"'' {q, e-t-iO) 
arises from the branch cut of f(x + iy) on the real axis for 
X > 1. Furthermore, using the fact that for |arg(— z)| < 



/(^) = ^(-)-'' 

sm(7r?7) 



rj z 

(3.45) 



we obtain in the regime — ti^g^j 3> 1 and ?/ <C 1 for 
the real part, 

Ref(^-^'^.'^)(g,e-f zO) = -{e^Hq) 



+ e{esq > ~3v\q\){e~vq)- 







+ 








e — vq 


e + vq 




Av 




2i 





+ Q{-3v\q\ > esq){e-vq) 



e — vq 



Av 



(3.46) 



In deriving this expression, we have neglected corrections 
of the order of to the prefactor. Note that the imagi- 
nary part given in Eq. (3.44) is proportional to ry, while no 
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such small prefactor appears in the real part, E g. ( 3^46). 
The first term on the right-hand side of Eq. ( 3.4E ) ex- 
actly cancels the inverse bare propagator ie — vq such 



that roo((7, e + iO) — (e — vq 
the dynamic scaling proper 



f(2- 



(g, e + zO) has 



r«,(q,e + *0) = |q|i-Vi(eM), 



(3.47) 



where the functions f^{x) and f^{x) refer to q > and 
q < respectively and {fi^{x))* = —f^{x). A graph of 
/+ (e/Cig) is shown in Fig. || for 77 — 0.2. 

Using the fact that for ^ we may identify 



Zq e 



Zq e 



-rjl 



(3.48) 



where the constant Zq = 0(1) has been absorbed by ^0, 
the corresponding spectral function can be written in the 
scaling form similar to Eq. (11. 



A{a{kF +p),u}) 



^0 



^oo(pC,Wt) 



(3.49) 



where r = S^/vp = 1/A, and the scaling function Aoo(9, e) 
is 



^00(9, e) 



Im 

TT 



1 



e-vq + ^ac 
1 + 



(g,e + «0) 







e — vq 


n+vq 





X < 





n/2 

+ 




V2" 


e-\~vq 


1 t — vq 


e — vq 




2 e+Dg 





for esg > u|(7| 

for v\q\ > eSq > —v\q\ 

for —v\q\ > eSq > — 3z;|g| 
for — 3-51(71 > eSq . 



(3.50) 



In the prefactor we have retained only the leading order 
in 77. For q = 0, the dynamic scaling function reduces to 



^oo(0,ej = -|e 



(3.51) 



If g ^ 0, the dynamic seeding function can be written in 
the dynamic scaling fornxj 



(3.52) 



where hoo{x) is defined via Eq. ( 3.5Q| ). A plot of the 
dynamic scaling function hooi^/vq) is given in Fig. I for 
r] = 0.2. For a comparison, we have also plotted the 
generally accepted scaling function of the TLM, given in 
Eq. (1.7). The qualitative agreement between the two 



plots becomes even better for smaller rj. Recall that Eq. 
( |3.50| ) is only vahd for j^l 1, |e| 1 and for large 
||e| — In this lim it th e small constant flao on the 

right-hand side of Eq. (3.50) can be neglected. In terms 
of the physical variables p — q/£, and uj — e/r — evp/^ 
the condition ||e| — v\q\\ ^ 1 becomes ||a;| — Vc\p\\ ^ 
1/^, where Vc — vvp. Using the fact that the function 
Aooiq, e) satisfies the scaling relation 



Aoo{sq, se) = s"^ ^^00(9, e) 



(3.53) 



we may write Aoo{p(,,LtJT) — t^~^ Aoc{vpp,uj). With 
= Aq'' we finally obtain for the spectral 
function at the Luttinger liquid fixed point. 



A{a{kp+p),u;) - A-"j\LU-vM''-' 







UJ — VcP 







X < 



4 
I 1 





rj/2 

+ 


1 LU — VcP 


r,/2" 


UJ — VcP 




2 uj-\-VcP 





for LOSp > Vc\p\ 

for Vc\p\ > ujSp > —Vc\p\ 

for -Vc\p\ > ujSp > -3vc\p\ 
for — 3wcb| > . 



(3.54) 
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^ 2 
I > 



- 






f Re/: 

im/: 



2 4 

E / yq 

FIG. 4: Real and imaginary part of /^(e/Cg) = r'cx)(q,e + 

(q,e + iO)l /Igl^-" for 7? = 



[(£ - wg) + f 



{2-fi,Z,u) , 



iO)/|g 

0.2, where f e + iO) is th e sub tracte d irr educible 

two-point scaling function, see Eqs. (|3.4^) and (3.46). 



Let us now compare Eq. ( 3.54 ) with the bosonization 
result for the TLM given i n Eg . (1.2). First of all, for 
p = we obtain from Eq. ( 3.54 ) 



(3.55) 



which corr espo nds to Eq. ( 3.52| ) and agrees exactly with 
the result (|l.2| ) obtained via bosonization. Thus, at least 
for p = 0, our truncation of the exact RG flow equations 
is reliable. Moreover, Eq. ( |3.54| ) exhibits the same scaling 
behavior Aj akp + asp, suj) = s^^^A{akp + ap, uj) as the 
expression ( |l.2| ) obtained via bosonization. For finite p, 
however, the detail ed li ne-shape of the spectral function 
predicted by Eq. (3.54) is different from Eq. (1.2). In 
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8 
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RG 
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e/ vq 



FIG. 5: Graph of the dynamic scaling func tion h og (e/ vq) = 
A^{q,e)/\q\''-^ for 77 = 0.2, see Eqs. (0|) and For 



smaller rj the RG result is indistinguishable from the bosoniza- 
tion result fiTLi^/vq) on this scale. 



particular, for lu 
predicts 



VcP + the bosonization result (|l.2| ) 



ATL{kF+p,u;) ^ Ap'^\2vcpr/^\u; ~ v,p\-^+^l^ 



while our RG calculation yields 

A{kF +p, cj) - Aq "27711^ - v^p\-^"^'^ 



(3.56) 



(3.57) 



which diverges with a different exponent than Eq. ( |3.56[) . 
Similarly, for uj — » ~VcP — the bosonization result (1.2) 
predicts 



ATLikF+P,u;) ^ A-''^|2t;,p|-i+''/2|c. -f v^''/^ 



while Eq. ( |3.54D yields 

A{kF + P,iu) ^ Ao'^2t]\2vcP\-^\u; - 



VcPl' 



(3.58) 



(3.59) 



Note, however, that Eq. ( |3.57 ) has been derived assum- 
ing — v^q^ I 1, which means in terms of the physical 
variables — Vc\p\\ ^ l/r = A. For A ^ this condi- 
tion is always satisfied unless precisely lu = iucp. This 
limiting procedure describes the spectral function pre- 
cisely at the critical point, because the length ^ = wf/A 
associated with the infrared cutoff can be sent to infinity 
before the limit uj — > VcP is taken. Obviously, there exists 
another regime where both |e| = \uj\t and |g| = \p\S, are 
large, but where 



v\q\ 



< 1 



(3.60) 



This corresponds to first taking the limit lu — )■ iVcP, 
and then removing the infrared cutoff A = 1/t — *■ 0. In 
this limit the spectral function behaves completely dif- 
ferent. Keeping in mind that in practice the TLM de- 
scribes a physical system only above some finite energy 
scale A* below which backscattering becomes important 
or a crossover to higher dimensionality sets in, this limit- 
ing procedure is physically relevant to describe quasi Id 
systems. If we naively u se ou r result for the irreducible 
self- energ y given in Eq. ( 3.41 ), we find with the help of 
Eq. (3.45) that for e ^ vq and for small 77 



-{ie — vq) . 



(3.61) 



Since by assumption \q\ ^ 1, this is a negligible correc- 
tion to the bare inverse propagator ie — vq. It follows 
that in this regime our scaling function is simply 



A^{q,e) w 6{e ~ vq) 



(3.62) 



so that the physical spectral function exhibits a small 
but still finite quasiparticle peak with a weight of the 
order of Zi^ ~ (A^/Aq)''. Note that this quasiparticle 
peak is due to the finite crossover energy scale A* and 
does not occur in the TLM at zero temperature where 
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the infr ared cutoff A can be reduced to zero. In deriv- 
ing Eq. ( 3.62| ) we have introduced the Dirac sea, ignored 
the non-Unearity in the energy dispersion as well as other 
irrelevant couplings, which will broaden the (5-function. 
In principle all these effects can be systematically taken 
into account within our formalism. The calculations are 
quite tedious and are beyond the s cope of this work. It 
seems, however, that in the regime (3.60) the behavior of 
the dynamic scaling function of a generic quasi Id inter- 
acting Fermi system with dominant forward scattering 
is not correctly described by the scaling function o f th e 
Tomonaga-Luttinger model A'Ti^{q,e) given in Eq. (1.6). 



IV. SUMMARY AND CONCLUSIONS 

In this work we have shown how to calculate 
the momentum- and frequency-dependent single-particle 
spectral function by mcaps of a particular version of the 
functional RG method.tJ For simplicity, we have applied 
this method to the exactly solvable TLM. Our trunca- 
tion scheme of the exact hierarchy of RG flow equations 
consists of the following steps: 

1. Calculate the flow of the relevant and marginal cou- 
plings perturbatively for small interactions and ad- 
just the initial values such that the RG flow ap- 
proaches a fixed point. 

2. Calculate the momentum- and frequency- 
dependent part of the four-point vertex (which 
involves an infinite number of irrelevant couplings) 
in powers of the marginal part of the four-point 
vertex. 

3. Calculate the anomalous dimension in powers of the 
marginal part of the four-point vertex. 

4. Substitute the results of steps 2. and 3. into the ex- 
act flow equation of the two-point vertex and then 
solve this equation exactly. 

The result for the spectral function A{k,u:) is quite 
encouraging: it has the correct scaling properties, and 
for k = ±fci? agrees with the exact bosonization result. 
For finite k — kp and w, we have found evidence that the 
spectral line shape in the vicinity of the Luttinger liquid 
fixed point exhibits some non-universal features. 

Our work also shows how non-Fermi liquid behavior 
in a strongly correlated Fermi system can be detected 
using modern functional RG methods. Recently several 
authors have presented numerical studies of the one-loop 
flow equations for the marginal part of the urwrescaled 
four-point vertex for Hubbard models in 2(i.LJ'LJ The 
marginal part is obtained by setting all frequencies equal 
to zero and projecting all Wiaye-vectors onto the Fermi 
surface. Ignoring ambiguitiesE3 related to different orders 
of limits, the unrescaled part of the four-point vertex is 



then approximated by 

rA^(ki, Wi,k2, W2;k2,t^2,ki, Wi) 



r\^\k'^i,0,k^2,0;kj.2,0,kj^i,0) . (4.1) 



,ext-nearest neighbor 
'^^ of the one-loop flow 



For the 2d Hubbard model w: 
hopping the numerical analysi; 
equations for the set of marginal couplings contained in 

the vertices rj(*''(k^;^, 0, k'^2i 0; 1^F2, 0, ki?i, 0) exhibits a 
runaway flow to strong coupling at a finite length scale 
U = ln(Ao/A,), which is typically of the order of 10^. 
The usual interpretation of this runaway flow is that it. 
signals some instability of the normal metallic stateE3'Ell 
In Ref. ^ we have pointed out that this runaway flow 
might also be an artifact of the one-loop approximation, 
and that the properly renormalized vertices could remain 
finite at the two-loop order. The crucial point is that the 
renormalized vertices involve wave-function renormaliza- 



tion factors (see Eq. (2.18)), which for a strongly corre- 
lated system can cancel a possibly strong enhancement 
of the unrescaled vertices. Our simple Id model allows 
us to study such a scenario in detail. In this case F^'' can 
be parameterized in terms of a single marg inal coupling 
gi, which is defined in analogy to Eqs. (2.6) and (2.41), 



T)^ je-i (a'l , 0, a2fc_F , 0; a2fc_F , 0, aifc_F , 0) 

— -^a[a'^;a2ai9l ■ (4-2) 

At the two-loop order, the flow equation of gi for our 
model is £3 



digi 



This implies 



91 = 



4 3i 



.90 



1 



(4.3) 



(4.4) 



Obviously, gi diverges at a finite scale U — 2/{vogo)^. 
However, this divergence and the associated runaway- 
fiow to strong coupling are unphysical, because the cou- 
pling gi in Eq. (4.2) is not the properly renormalized cou- 
pling that can be identified with the usual (72-interaction 
of the TLM. The latter is defined as a model describ- 
ing the fixed point of the RG. Thus, the g2-coupling that 
appears in the TLM should be identified with the renor- 
malized coupling at the RG fixed point, which is related 
to the fixed point value of our rescaled coupling gi, 



vog2 = ^fim gi = Yaa [zfi/ogi] 



(4.5) 



see Eqs. ( 2.18 ) and ( 2.41 ). Similar relations between ver- 
tex functions and interaction parameters at the RG fixed, 
point are well known from Landau Fermi liquid theory.Ej 
In our simple model with a linearized energy dispersion 
the two- loop flow equation of the rescaled coupling is sim- 
ply digi = 0, so that the limit / — > 00 in Eq. (4.5) indeed 
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exists. Note that the vanishing wave-function renormal- 
ization factor Zi exactly compensates the diverging un- 
rescaled couphng gi such that the rescaled coupling gi 
remains finite and small. We believe that the above inter- 
pretation of the runaway flow of the RG-flow equations 
for the marginal part of the unrescaled vertices V^f^ is 
not only specific to Id, where the RG fixed point does 
not correspond to a Fermi liquid. Assuming that in 2d 
the RG fixed point corresponds to a strongly correlated 
Fermi liquid with <g; 1 , we expect that the unrescaled 
vertices given in Eq. (O) will flow to a finite but large 
value, which numerically might be indistinguishable from 
a runaway flow to infinity. At the same time, the prop- 
erly rescaled vertices f [''^ defined in Eq. (2.18) can remain 
finite and small. 

The problem of calculating correlation functions of 
many-body systems at finite wave-vectors or frequen- 
cies using RG methods has not received much attention. 
For classical systems, functional RG calculations of the 
momentum-dependence of correlation functions can be. 
found in the textbook by Ivanchenko and Lisyansky.E3 



Here we have presented a functional RG calculation of 
a momentum- and frequency-dependent correlation func- 
tion of a non-trivial quantum-mechanical many-body sys- 
tem. For a special 2d system a similar calculation has 
recently been performed by Ferraz,t£l who used the field 
theory RG. We believe that the method described in this 
work will also be useful to study other problems where 
no exact solutions are available. For example, the spec- 
tral function of one-dimensional Fermi systems where 
backscattering or Umklapp scattering are relevant cannot 
be calculated exactly by means of bosonization or other 
methods. Using our RG method, it should be possible 
to obtain the spectral function of Luther-Emery liquids 
even away from the Luther-Emery point. 
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